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Constraining scalar dark matter with Big Bang nucleosynthesis and atomic 

spectroscopy 

Y. V. Stadnilfl and V. V. Flambaunfl 

School of Physics, University of New South Wales, Sydney 2052, Australia 
(Dated: April 14, 2015) 

Scalar dark matter can interact with Standard Model (SM) particles, altering the fundamental 
constants of Nature in the process. Changes in the fundamental constants during and prior to Big 
Bang nucleosynthesis (BBN) produce changes in the primordial abundances of the light elements. 

By comparing the measured and calculated (within the SM) primordial abundance of ^He, which 
is predominantly determined by the ratio of the neutron-proton mass difference to freeze-out tem¬ 
perature at the time of weak interaction freeze-out prior to BBN, we are able to derive stringent 
constraints on the mass of a scalar dark matter particle cj) together with its interactions with the 
photon, light quarks and massive vector bosons via quadratic couplings in 0, as well as its inter¬ 
actions with massive vector bosons via linear couplings in 0. We also derive a stringent constraint 
on the quadratic interaction of (j) with the photon from recent atomic dysprosium spectroscopy 
measurements. 


I. INTRODUCTION 

Measurements m and Standard Model (SM) cal¬ 
culations Hi] of the primordial abundances of the 
light elements produced during Big Bang nucleosynthe¬ 
sis (BBN) provide a valuable probe into physics beyond 
the SM, including axion-like pseudoscalar dark matter 
(DM) [To], neutrinos and relativistic non-SM particles 
in the early Universe mini, DM reactions during BBN 
m na nsj, variations of the fundamental constants of 
Nature m 0 0 0 13 ehui and supersymmetric mod¬ 
els [22]. BBN considerations have also been invoked to 
explain the fine-tuning of the weak scale in the SM [23] . 
For contemporary reviews of various aspects of BBN, we 
refer the reader to Refs. [2411^ . 

In our previous work [26] , we pointed out that a global 
cosmological evolution of the fundamental constants due 
to an oscillating scalar DM field (f that interacts with SM 
particles via quadratic couplings in (f is most stringently 
constrained by determination of the ratio of the neutron- 
proton mass difference to freeze-out temperature at the 
time of the weak interaction freeze-out prior to BBN. 
The stringency of the constraints arises due to the high 
energy density of DM in the early universe compared 
with that at present, which enhances the effects of any 
possible variations of the fundamental constants due to 
the underlying field (j). 

In the present letter, we derive constraints on the mass 
of a scalar DM particle (j) together with its interactions 
with the photon, light quarks and massive vector bosons 
via quadratic couplings in 0, as well as its interactions 
with massive vector bosons via linear couplings in 0, by 
comparing the measured and calculated (within the SM) 
primordial abundance of ^He, which is predominantly de- 
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termined by the ratio of the neutron-proton mass differ¬ 
ence to freeze-out temperature at the time of weak in¬ 
teraction freeze-out prior to BBN. We also derive a con¬ 
straint on the quadratic interaction of (f with the pho¬ 
ton from recent atomic dysprosium spectroscopy mea¬ 
surements. Our constraints are found to rule out large 
regions of previously unconstrained scalar DM parameter 
space. 


II. THEORY 

In the present work, we consider a scalar DM field 0, 
which can couple to the SM fields via the following linear- 
in-0 interactions: 


= ( 1 ) 

/ ^ 

where the sum runs over all SM fermions /, rrif is the 
standard mass of the fermion, / is the fermion Dirac field 
and f = p'y°, 
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1 

int 


(!> 
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( 2 ) 


where are the components of the electromagnetic 
field tensor, and 
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int 




h.y 2 
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(3) 


where the sum runs over all SM massive vector bosons 
U, My is the standard mass of the boson and are 
the components of the wavefunction of the correspond¬ 
ing massive vector boson. Ax is a very large energy 
scale, which is strongly constrained by equivalence prin¬ 
ciple tests, including lunar laser ranging [27] [28] and the 
EotWash experiment [29] [30] (see also [31] for constraints 
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from stellar energy loss bounds). Eqs. 0 , 0 and 0 
alter the fundamental constants as follows, respectively: 


rrif rrif 



( 4 ) 


1 =F 




( 5 ) 


Mv Mv 



(6) 


The scalar DM field (j) can also couple to the SM fields 
via quadratic-in -0 interactions, with the replacement 
(fy/kx {(k/k'xf in Eqs. 0-0. is a large energy 
scale, which is constrained by astrophysical observations, 
most notably bounds from supernova energy loss, and 
equivalence principle tests [3^ . 

Scalar DM can thus induce a cosmological evolution of 
the fundamental constants. Changes in the fundamental 
constants during and prior to BBN alter the primordial 
abundances of the light elements. There are two limiting 
cases of particular interest to consider: 

(i) 0 is an oscillating field and interacts via SM par¬ 
ticles via quadratic couplings in 0. This occurs when 
rricf) where is the mass of the scalar DM par¬ 

ticle and H{t) = l/2t is the Hubble parameter as a func¬ 
tion of time [33] . The energy density of a non-relativistic 
oscillating scalar DM field is given by Pscaiar — (0^) 

and for a non-relativistic cold field evolves according to 
the relation 


III. RESULTS 

The loosest constraints on the parameters of 0 are 
obtained by assuming that 0 is a non-relativistic cold 
field at all times, which is produced either non-thermally 
(through vacuum decay [35]) or thermally (with a very 
large mass). An example of a cold scalar DM field is a 
classical oscillating condensate, 0 = 0o cos{ujt), which os¬ 
cillates with frequency uj ~ rriffyC^/h. We note that BBN 
and Cosmic Microwave Background (CMB) measure¬ 
ments do not rule out the existence of a relativistic scalar 
DM field [1] , and that a relativistic scalar DM field should 
provide even more stringent constraints on the underly¬ 
ing parameters than those derived in the present work, 
since for a relativistic oscillating DM field, the mean DM 
energy density evolves according to pdm oc [1 z{t)]^. 

Constraints on scalar DM parameters follow from mea¬ 
surements and SM calculations of the abundance of ^He 
produced during BBN, which depends strongly on the 
neutron-to-proton ratio at the time of the weak interac¬ 
tion freeze-out (Tp = {Ow)/^ 

0.75 MeV [33]), where Tp is the weak interaction freeze- 
out temperature, Ow is the Weinberg angle, a is the 
electromagnetic fine-structure constant, Mpianck is the 
Planck mass and 6 is a numerical constant. The neutron- 
to-proton ratio at the time of weak interaction freeze-out 
is given by 


n 

P 


— g QupIT'f 

weak 


( 9 ) 


Pbm = 1-3 X 10 ^ [1 -h ( 7 ) 

where z{t) is the redshift parameter and the present mean 
DM energy density is determined from WMAP measure¬ 
ments [34]. 

(ii) 0 is a non-oscillating field and interacts with SM 
particles via linear or quadratic couplings in 0. This 
occurs when ^ H{t). The energy density of a 
non-oscillating scalar DM field is given by Pscaiar = 
(0^) /2 and, due to Hubble friction, is approximately 
constant while the field remains non-oscillating: 

pi AT 

PDM ~ 1-3 X 10 “^ [1 + z{tm)f - J, ( 8 ) 

where z{tm) is defined by H{tm) = '^ 0 - 

For constraining the parameters of scalar DM from 
measurements and SM calculations of the primordial 
abundance of ^He, it suffices to consider only the effects 
of variation of the fundamental constants due to scalar 
DM on the ratio of the neutron-proton mass difference 
to freeze-out temperature at the time of weak interac¬ 
tion freeze-out (tp ~ 1.1 s), which determines the abun¬ 
dance of neutrons available for BBN (the vast majority 
of these neutrons are ultimately locked up in ^He). The 
corresponding range of scalar DM particle masses for (i) 
is hence ^ 10 “^^ eV, while the corresponding range 
of scalar DM particle masses for (ii) is ^ 10“^^ eV. 


where Qnp is the neutron-proton mass difference: 

Qnp = mn-mp = aaAqcD + {rrid - m^), ( 10 ) 

with the present-day values (aQfAQCD)o = —0.76 MeV, 
where Aqcd is the Quantum Chromodynamics scale and 
a is a numerical constant, and {rrid — mu)o = 2.05 MeV 
[36] . From the measured and predicted (within the SM) 
primordial ^He abundance, 10®"^P(^He) = 0.2474 ± 0.0028 
[ 2 | and (^He) = 0.2486 ± 0.0002 [ 8 |, we find the 

measured and predicted n/p ratio at the time of BBN 
freeze-out to be 


exp 

= 0.1420 ±0.0021, (11) 

BBN 

theor 

= 0.1428 ±0.0001. (12) 

BBN 

Extrapolating back to the time of weak interaction freeze- 
out (At 180 s) with a neutron half-life of = 880 s 
[34] . gives the measured and predicted n/p ratio at the 
time of weak interaction freeze-out: 

exp 

= 0.1801 ±0.0026, (13) 

weak 
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n 

P 


theor 

= 0.1811 ±0.0002, 

weak 


where we have made use of the relation cos(6>w) = 
^ M\y/M z — 0.882 [34]. 


which correspond to the following measured and pre¬ 
dicted Qnpt'^F ratio: 

1,714 ±0.015, (15) 


The constraint (18) can be expressed in terms of (j) 
and Ax as follows (retaining only variations in the fun¬ 
damental constants that are induced by the linear-in-0 
interactions considered in the present work): 


I Qnp A 
\Tf / 


1.709 ±0.001. 


(16) 


The relative difference in the measured and predicted 
QupITf ratio is found to be (adding the experimental 
and theoretical uncertainties in quadrature): 

A{Qnp/TF) ^ Q_gQ 33 ^ 0.0085. (17) 

Qnp/T'p 

can be interpretted as a constraint on temporal vari¬ 
ations in the underlying fundamental constants from the 
time of weak interaction freeze-out until the present time: 
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= 0.0033 ± 0.0085, 
(19) 


0.08- 
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1.59 
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A(md - rriu) 

{rrid - rriu) Mw 

AAqcd ^ 1 AMpianck 
AqcD 3 Mpianck 


3,32±0r_4,65^ 
Mz 


= 0.0033 ± 0.0085, 
(18) 


where kx = ±1 correspond to the relevant signs in the 
Lagrangians 0-(i- For a non-oscillating scalar DM 
field (rricf) 10“^^ eV at tp ~ 1.1 s), for which the energy 
density evolves according to (|^, this leads to 


J_| 


3/4 

"0.08/^^ 1 1.59 

f i^drrid 

KuTTlu^ 

^ 3.32k\y 


m0 

iv3 X 10-16 eV, 

) 

_ A^ rrid- rriu 

y Ad 

Au J 

± Aw 

Az 


(0.4 ± 1.0) X 10“^^ eV“^ 

( 20 ) 


where we have made use of the fact that | (0)weak I ^ 

I (0)q I and the relation [1 ± z{tm)]/{l ± ^f) - 
during and after BBN (but at much earlier times than 
electron-proton recombination), and assumed that scalar 
DM saturates the present-day DM energy density. We 
note that a single type of measurement does not give con¬ 
straints on the individual parameters appearing in (20), 
but rather gives constraints on a combination thereof. 
However, we can extract useful information about the 
underlying parameters by sequentially assuming that in¬ 
dividual terms within (20) dominate the others. The 
resulting region of parameter space excluded by compar¬ 
ison of measurements and SM calculations of the primor¬ 
dial abundance of ^He is shown in Fig. [T] in terms of 
the parameter Ay = |AvfA^/(A^ — 1.40Avf)|, assuming 
Kw = For the linear interaction of 0 with the photon 
and light quarks, the BBN limits are weaker than exist¬ 
ing limits from equivalence principle tests Ei EO] and 
dysprosium spectroscopy measurements [37] in the mass 


range of interest. 

The constraint ( [l8| ) can also be expressed in terms of 
0 and A(j^ as follows (retaining only variations in the 
fundamental constants that are induced by the quadratic- 
in-0 interactions considered in the present work): 


)weak 




ro.osy 1.59 

- - H- 

_ (A !;^)2 TUd - rriu 

3.32K'jy Amn'z' 


( K'^rtid <m„ \ 

bA')2 (AOV 
= 0.0033 ± 0.0085. 
( 21 ) 


For an oscillating scalar DM field {m^ ^ 10 eV at 
tp ~ 1.1 s), for which the energy density evolves accord- 
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ing to Q, this leads to 


1 




'0.08/^^y 

.7^ 


1.59 / n'^nud i^'yjnvX 3.32/^'^ 

md-mu Ua')2 (A'^)2 

-^] ==(1.0±2.5)xlO-“eV-‘, 

( 22 ) 


FIG. 1. (Color online) Region of scalar dark matter param¬ 
eter space ruled out for the linear interaction of 0 with mas¬ 
sive vector bosons. Region below blue line corresponds to 
constraints derived in the present work from comparison of 
measurements and SM calculations of the ratio Qnp/Fp, for 
the case Ky ^ A^, Ky < |A^iAc^(md-r7^n)/(A^^?nc^ - A^mn)!. 


where we have made use of the fact that ^ 

and assumed that scalar DM saturates the present- 
day DM energy density. For a non-oscillating scalar DM 
field (rriff) <C 10“^^ eV at tp ~ 1.1 s), for which the energy 
density evolves according to (|^, this leads to 


J_l 

( rn^ ^ 

3/2 

'0.08k'^ 1.59 



\ 1 3.32/1:^ 

4.65k'z' 


(3 X 10-16 eV, 

) 

_ (A^)^ rud - m„ 



' (AW)^ 

to 


(23) 


where we have made use of the fact that ^ 

and the relation [1 + z{tm)]/{l + ^f) - V^fT^ 
during and after BBN (but at much earlier times 
than electron-proton recombination), and assumed that 
scalar DM saturates the present-day DM energy den¬ 
sity. The resulting regions of parameter space ex¬ 
cluded by comparison of measurements and SM cal¬ 
culations of the primordial abundance of ^He are 
shown in Figs. 0-0 in terms of the parameters 
(A;)2 = |(A;)2(A')2(m,-m„)/[(A;)2m,-(A')2m„]| 
and (AV)2 = |(A'^)2(A'^)V[(A'^)2-1.40(A'^^)2]|, assum- 
ing and Note that we have extrapo¬ 

lated the constraints from the two limiting mass ranges 
to the connecting central region. 


We have also derived constraints on the quadratic in¬ 
teraction of <p with the photon, using data from re¬ 
cent atomic dysprosium spectroscopy measurements in 
Ref. [37] pertaining to oscillating variations in a. These 
constraints are presented in Fig. 


IV. CONCLUSIONS 

By comparing the measured and calculated (within the 
Standard Model) primordial abundance of ^He, which is 
predominantly determined by the ratio of the neutron- 
proton mass difference to freeze-out temperature at the 
time of weak interaction freeze-out prior to Big Bang nu¬ 
cleosynthesis, we have derived stringent constraints on 
the mass of a scalar dark matter particle (p together with 
its interactions with the photon, light quarks and mas¬ 
sive vector bosons via quadratic couplings in 0, as well 
as its interactions with massive vector bosons via lin¬ 
ear couplings in p. We have also derived constraints 
on the quadratic interaction of p with the photon from 
recent atomic dysprosium spectroscopy measurements. 
The new constraints derived in the present work rule 
out large regions of previously unconstrained scalar dark 
matter parameter space. 
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FIG. 2. (Color online) Region of scalar dark matter parameter 
space ruled out for the quadratic interaction of with the pho¬ 
ton. Region below blue line corresponds to constraints derived 
in the present work from comparison of measurements and 
SM calculations of the ratio Qnp/Fp, for the case Ag, 

^ A'^. Region below black line corresponds to constraints 
from supernova energy loss bounds [32] . Region below red line 
corresponds to constraints derived in the present work from 
the data of recent atomic dysprosium spectroscopy measure¬ 
ments in I3Z1- 



FIG. 3. (Color online) Region of scalar dark matter parame¬ 
ter space ruled out for the quadratic interaction of 0 with 
light quarks. Region below blue line corresponds to con¬ 
straints derived in the present work from comparison of mea¬ 
surements and SM calculations of the ratio Qnp/Fp, for the 
case Aq Aq Ay. Region below dashed black line 

corresponds to constraints on A^ from supernova energy loss 
bounds [32]. Strictly, the limits are independent, since they 
contain different linear combinations of the quark interaction 
parameters. 
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